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Abstract: We systematically analyze the integrability of a Pauli system in Lorentz violat- 
ing background at the non-relativistic level both in two- and three-dimensions. We consider 
the non-relativistic limit of the Dirac equation from the QED sector of the so-called Stan- 
dard Model Extension by keeping only two types of background couplings, the vector 
and the axial vector b^. We show that the spin-orbit interaction comes as a higher order 
correction in the non-relativistic limit of the Dirac equation. Such an interaction allows the 
inclusion of spin degree non-trivially, and if Lorentz violating terms are allowed, they might 
be comparable under special circumstances. By including all possible first-order derivative 
terms and considering the cases a 7^ 0, b 7^ 0, and bo 7= one at a time, we determine the 
possible forms of constants of motion operator, and discuss the existence or continuity of 
integrability due to Lorentz violating background. 
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1. Introduction 

In classical mechanics, integrability of a Hamiltonian system can be defined as the exis- 
tence of a set of n functionally independent constants of motion X±, X2, • X n including 
the Hamiltonian itself if the system has n degrees of freedom. These constants of mo- 
tion must be in involution. The concept of integrability can be extended to so-called 
superintegrability which requires at least one additional constant of motion (minimally 
superintegrable) and allows totally n — 1 additional constants of motion (maximally su- 
perintegrable) 1 . These additional constants of motion, however, are not necessarily in 
involution among each other, nor with X±,X2, In a very similar manner, the no- 

tion of (super) integrability can also be defined in quantum mechanics through well-defined 
linear constants of motion operators which are now supposed to be algebraically indepen- 
dent [1, 2, 3, 4]. In quantum mechanics, integrability not only simplifies the calculation of 
energy levels and wave functions, but also provides a complete set of quantum numbers, 
which characterize the system completely. Superintegrability, on the other hand, may en- 
tail exact solvability. The harmonic oscillator [2], and the Kepler, or Coulomb system [1] 
are well known examples of superintegrability. 

1 This is due to the fact that the phase space is 2n dimensional and at least one free degree of freedom 
has to be left in order to have dynamics in the system. 
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Systematic searching for superintegrability in Hamiltonian systems with kinetic en- 
ergy quadratic in momenta and with only coordinate dependent potentials was started 
quite some time ago [3]. In these works the constants of motion were considered up to 
second-order in momenta in order to search for quadratic (super)integrability Both in 
classical and quantum mechanics, second-order integrability of such Hamiltonians mim- 
ics the separations of variables which is the first step for exact solvability. Such close 
connection doesn't hold for the systems with velocity dependent potentials V(p,A) [5]. 

Recently, the (super)integrability of spin-dependent Hamiltonian systems with a generic 
scalar potential has been studied in a systematic way up to first-order [6]. There are other 
known systems which introduce velocity-dependent interactions without requiring an ex- 
ternal vector potential A. For example, if one allows violation of Lorentz symmetry at the 
relativistic level, the extra fields which describe such violation can also induce velocity- 
dependent interactions at non-relativistic level. The aim of this study is to systemati- 
cally analyze (super) integrability of a Hamiltonian system (more precisely the hydrogen 
atom) with spin-orbit interactions in some Lorentz violating backgrounds both in 2- and 
3-dimensions. 

There exists a well established framework [7], called Standard- Model Extension (SME), 
for the study of the Lorentz and Charge-Parity-Time reversal (CPT) violation. The Stan- 
dard Model Extension is the generalization of the Standard Model (SM) with additional 
Lorentz and CPT violating interactions introduced through some tensorial background 
fields. It can be considered a low-energy limit of a fundamental theory in which the 
Lorentz and CPT symmetries are exact but broken spontaneously when evolved down to 
low energy scale, due to the existence of these tensorial background fields. SME is one of 
the elegant way to formulate the problem. Of course the coordinate independence of phys- 
ical observables are maintained by requiring observer Lorentz invariance which describes 
the transformations of coordinates. Lorentz violation considered in the model is so called 
particle Lorentz violation which describes rotations and boosts of particles and localized 
fields but not background fields in a specific observer's inertial frame. Formulating the 
model in non-Minkowski spacetimes [8] leads to spacetime dependent coupling coefficients. 
Our approach here is minimalistic, i.e., constant coefficients. 

In this study we consider only the effects of vectorial background couplings and 6 M to 
the (super)integrability of what we call the Pauli system. Here a M and 6 M are related to the 
vacuum expectation values of some vectorial background fields. Seeking (super)integrability 
for a Pauli system in such LV background and aiming exact solvability could be of some 
interest under certain circumstances. For example, the current bound on &o [11], the time 
component of b^, is 6q 1 > 1CP 3 cm and it is much weaker than that on |6|, the space part 
of the bp. So, for a system with timelike 6 M coupling the perturbative approach fails if the 
effective size of the system is of the order of 10~ 3 cm, and then the search for an exact 
solution becomes unavoidable. 

The outline of the paper is as follows. In the next section, Section 2, we introduce 
the relevant part of the SME Lagrangian and consider the non-relativistic limit. Then, in 
Section 3, the 2-dimensional Pauli system in LV background is discussed in three separate 
cases. In Section 4, before discussing the generalization of the problem to 3-dimensions, 



- 2 - 



we present the determining equations in the most general form. We then shortly discuss 
the perturbative approaches to the system in Section 5 and conclude in Section 6. 



2. The model and its non-relativistic limit 

There exist vast number of studies on the SME both at theoretical and phenomenological 
levels [9, 10]. Here we work in the QED sector of the SME for electron [12, 13, 14]. The 
QED Lagrangian for electron can be given as 

^electron = ^ ^ B M V ~ , 

r M = 7^ + ca m 7 A + d Xfl j57 X + ... 

M = m + arf + b^-f + \h^ u , (2.1) 

where = + ieA^ with the vector potential (/ D M g = fD^g - {D^f)g), a M 
and 6 M are the (pseudo-)vectorial CPT-odd coefficients and the other terms are CPT-even 
tensorial coefficients. These are the only terms that are obtained from the SME but 
can have more terms originating from non-renormalizable higher-dimensional operators. 

It is clear from the above Lagrangian that LV doesn't require CPT violation but the 
vice versa is true [15]. In this study, we consider the only two CPT-odd vectorial couplings 
and b^ and neglect the tensorial ones 2 . A comprehensive analysis of field redefinitions 
and redundant parameters of the model can be found in Ref. [16]. For example, the 
term can be absorbed by redefining the fermion field by a phase factor exp{ia^x^) and 
thus the spectrum will be unaffected by such term. Depending on the complexity of the 
model (photon interaction, fermion mixings, etc.), this situation could change. A similar 
transformation can be found for the case if, for example, a free massless fermion is 
assumed, but in general the b^ term is non-trivial. So, in our discussion we keep both a M 
and b^ terms. 

We now consider the non-relativistic limit of the Dirac equation for electron starting 
from Eq. (2.1) with LV vectorial couplings only. We take the external electromagnetic field 
as = (A , 0). The Dirac equation becomes 3 

{Puff - eA^ - - - m)^ = , (2.2) 

where we adopt the Dirac representation for the 7 matrices. For convenience one can define 
the energy of the electron po = p' + m (by extracting the rest energy) so that p' becomes 
much smaller than m in the non-relativistic limit. If we also assume eAo, \a\, \b\ <C m, we 
obtain the following equation for the 2x1 spinor ?/>i, which is the large component of 1/}, 



{p' -m-aQ-er-b)- —(cr ■ (p - a) - b ) A 1 (er • (p - a) - b ) 



y>i = o, 



A=[l + P °- eA °^-"- b ), (2.3) 



2 Note that allowing non-renormalizable higher-dimensional operators induces additional vectorial fields. 
3 We use the natural units, h = 1, c = 1, throughout the paper. 
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where <x = (a±, (T2, 03) are the Pauli matrices. Keeping the first two terms in A 1 and 
making some straightforward manipulations lead to the following equation 

d 
at 

1 4 

D 6 it %£ 

Hnr = — + eA - —r a Exp- + — — E ■ p + H LW , 

2m 4m z om J 4m z 

i?LV = -— ap + cr-6 -cr p + p 2 + - — ^ cr per b, (2.4) 

m m 4m^ 4m z 

where E is the electric field and we keep up to linear terms in LV parameters, since they are 
presumed small, and we also neglect the constant terms. In i?NR> the first order correction 
terms to the usual Hamiltonian are the spin-orbit, relativistic, and the potential energy 
correction, respectively. We assume that the potential energy is spherically symmetric, 
Vb(r) = eAo(r) so that the term a ■ E x p with E = leads to the usual spin-orbit 

interaction cr • L with L the angular momentum operator. 

Had we neglected all these three corrections and considered the Lorentz symmetry 
preserved, there would be no spin dependence in the system (since we set the external 
magnetic field zero) and the spin degree of the electron would become trivial. The only 
spin dependence is through <r • L and we intend to keep that term and neglect the other 
two whose presence would considerably complicate our study. As seen from Eq. (2.4), Hjy 
also brings spin dependent interactions only through b^. The last two terms of Hvj come 
from the second term in A -1 expansion and are higher order with respect to the first three 
terms of i?LV- Hence, we neglect them in the rest of our calculation. We further set m = 1 
for simplicity. 

Taking into account all of the above remarks, and inserting the Coulomb potential of 
the nucleus, the Lorentz violating Hamiltonian -£/nr for the hydrogen atom follows from 
Eq. (2.4) as 

H = - {pl+pl+pfj +V (r) + V 1 {r)cr-L-a-p + cr-b-b cr-p, (2.5) 

where Vq = Ti/r, V\ = F^/V 3 with Y\ = — e 2 and T2 = — e 2 /4. In the absence of the 
LV terms we call the above Hamiltonian the Pauli system, even though a system with an 
external magnetic field is usually considered as Pauli system. The properties of the LV 
terms under the Charge (C), Parity (P), Time-reversal (T), and rotational symmetries are 
given in Table 1 

Table 1: Properties of LV terms in the Hamiltonian Eq. (2.5) under C, P, T, and rotational 
symmetries. 
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In order to investigate the first-order integrability and/or superintegrability we analyze 
the commutativity of a general first-order integral of motion 

X = (A + A ■ a)pi + (B + B- <t) P 2 + (C + C ■ a)p 3 + $ + * . a 

- % - {(A + A ■ a) x + (B + B ■ a) y + (Co + C ■ cr) z } , (2.6) 

with the Hamiltonian given in Eq. (2.5). In Eq. (2.6), A, B, C and <& where A = 
(Ai, A 2 , A 3 ) (B, C and <& are defined similarly) are real functions of r. The commutativity 
relation 

[H,X]=0, (2.7) 

has second-, first-, and zeroth-order terms in the momenta. Setting the coefficients of each 
of these terms equal to zero we get equations determining Aq, Bq, Co, <1?o and Ai, Bi, Cj, 
$i (i = 1,2,3). 

3. The analysis of 2D Pauli system in LV background 

In this section we consider the integrability of 2D Pauli system in LV background. There 
are three LV terms in the Hamiltonian given in Eq. (2.5). For the sake of simplicity we 
consider the effects of each term one at a time; the case with non-zero vector coupling 
(ai 7^ 0), with pure spacelike axial- vector coupling (bi ^ 0), and finally with pure timelike 
axial- vector coupling (b ^ 0). 

3.1 Lorentz violation with vector coupling (ai ^ 0) : 

In the absence of the axial-vector coupling terms in the Hamiltonian, Eq. (2.5) in two- 
dimensions is given as 

h zd ^ = \(pI +Pt) + Vo(p) + Vi(p) ^3 - a lPl - a 2P2 , (3.1) 

where again Vo = T\/ p and V\ = T2/p 3 , T\ and T 2 are constants defined in the previous 
section. Since now the Hamiltonian Eq. (3.1) is a diagonal matrix we can also choose X 
given in Eq. (2.6) as diagonal too and write 

X^ 0) = (Aq + A 3 a 3 )pi + (B + £3^73 )p 2 + $ + $3<X3 

~ {(A + A 3 a 3 ) x + (B + B 3 a 3 ) y } . (3.2) 

The requirement of vanishing of the commutator [H^^°\ X^^ ^] gives us a total of 12 
equations for A v , B u and (v = 0,3). The technique is to start from the highest order 
terms (second-order in our problem) and determine some functions in X^^ and then 
use these solutions to apply the same procedure repeatedly for the remaining orders until 
all unknown functions in X^^ are fixed. So, let's start from the second-order. We get 
the following six partial differential equations 4 

A vx = 0, B uy = 0, A vy + B ux = 0, (i/ = 0,3), (3.3) 

4 Throughout the paper, the subscripts x, y and z represent the partial derivatives with respect to the 
cartesian coordinates. 
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which can immediately be integrated to give 

A u = £ v + w v y , B v = rj v - w u x . (3.4) 

After introducing Eq. (3.4) into the coefficients of the first-order terms and separating 
the real and imaginary parts we have 4 determining equations for & u 

®ux = S v>3 - K [(w K x - rj K )yVi y - (w K y + £, K )yV lx - r) K V{\ + a 2 w v , 

= 5 v ,s-K[{w K y + € K )xV lx - (w K x - rj K )xV ly + ^Vi] - a\w v , {v = 0, 3) , (3.5) 

where 5 U>K is the Kronecker delta function. In order to satisfy the compatibility conditions 
for $j, we must have = and rj v = 0, v = 0, 3. Then, Eq. (3.5) can be easily integrated 
to give 

$ u = w u (a 2 x - aiy) , {y = 0, 3) , (3.6) 

where we used V\ = T 2 /p 3 and set the integration constants equal to zero without loss of 
generality. 

Two of the 12 determining equations, which are indeed the coefficients of the zeroth- 
order terms, have not been used up to now. If we introduce all the information gathered 
from the coefficients of the higher-order terms into the coefficients of the zeroth-order terms 
we get the following two conditions to be satisfied 

(aix + a,2y)V\w v = , (^ = 0,3). (3.7) 

Clearly we must have either (1) both a\ = and a 2 = 0, the components of the LV vector 
a M on the plane or (2) w v = 0. 

Case (1) corresponds to a LV vector of the form = (00,0,0,03) or even a pure 
timclikc vector (03 = 0). The latter is trivial since the system still remains Pauli system 
up to some constants. The former is non-trivial since a non-zero space component 03 
is allowed but it doesn't couple with the electron in xy plane. Therefore, under these 
conditions only, the symmetries of the Pauli system are restored. For example, we get the 
constant of motion as 

X < ^^°\a 1 = = a 2 ) = -(w + w 3 a 3 )(xp 2 - ypi) 

= -(w + w 3 a 3 )L 3 , (3.8) 

where u>o and w 3 are arbitrary constants. As expected both L 3 <gil and L 3 (g>a 3 are constants 
of motion. However, they are not different from each other since o 3 trivially commutes 
with the Hamiltonian in 2D. Therefore, there is only one extra constant of motion other 
than the Hamiltonian and there exists first-order integrability. 

For the Case (2), which is indeed the non-trivial one, X { 2 ^ 0) vanishes as seen from 
Eq. (3.8) and there is no first-order integrability. 
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3.2 Lorentz violation with pure spacelike axial- vector coupling (bi ^ 0): 

If we only have pure spacelike axial- vector coupling, i.e., bi 7^ 0, then in two-dimensions 
the Hamiltonian Eq. (2.5) and the integral of motion Eq. (2.6) are given as 

H m^ 0) = \(p\ +Pt) + MP) + Vi(p)<7 3 L 3 + *-b, (3.9) 
X^ 0) = (A + A ■ tr)pi + (5 + # • <r)p 2 + $ + * • o- 

- l -{(A + A-a) x + (B + B-a) y }, (3.10) 

where Vb = Ti/p and Vi = T 2 /p 3 are assumed. Since the Hamiltonian Eq. (3.9) is not 
diagonal, we can no longer choose the integral of motion operator x!jf£^ diagonal. 

In this case, unlike the previous case, there are 12 determining equations coming from 
the coefficients of the second-order terms since we have non-diagonal entries which double 
the number of equations. The ones from the diagonal elements are exactly the same as Eq. 
(3.3) and hence the solutions for A v and B u , (y = 0, 3) as given in Eq. (3.4), are the same. 
From the off-diagonal elements we also have the following 6 equations 

2yA 1 V 1 + A 2x = 0, 2yA 2 V l - A lx = , 
2x52V"! + B ly = , 2xBiVl - B 2y = , 
2V 1 (xA 2 - yB 2 ) + A ly + B lx = , 2V 1 (xA 1 - y B x ) - A 2y + B 2x = . (3.11) 

However, there is only a so called trivial solution (keeping in mind that V\ = T 2 / p 3 ) for 
the system of equations given in Eq. (3.11), which is A\ = A 2 = B\ = B 2 = 0. This has 
been shown both by hand and by using the Maple software [17]. 

When we introduce all of this information into the rest of the determining equations 
we get, in addition to the equations given in Eq. (3.5) (with of course a\ = and a 2 = 0), 
the following four equations for the contribution of the coefficients of the first-order terms 

2b 2 w 3 y + 2yVi$ 2 - $i* = , 2b 1 w 3 y + 2yV 1 $ 1 + <S> 2x = , 

2b 2 w 3 x + 2xV^ 2 + <& ly = , 26i w 3 x + 2zVi$i - <S> 2y = . (3.12) 

These equations can be solved and we find that £ v = 0, r\ v = 0, <I>o = (y = 0, 3) are 
solutions if either of the following conditions is satisfied 

• Case (1): b ± = = b 2 

• Case (2): w 3 = . 

Note that the determining equations coming from the coefficients of the zeroth-order terms 
are identically satisfied in either case. The constant of motion operator is the same as given 
in Eq. (3.6) but, since this time there is no condition on wo, L 3 ® I becomes a constant of 
motion. It is also important to note here that since the Hamiltonian Eq. (3.9) has a\ and 
a 2 in it, the fact that o~ 3 commutes with the Hamiltonian is no longer trivial. Hence, the 
term L 3 ® a 3 that comes with w 3 should now be considered as independent. Obviously, in 
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the presence of b\ and b 2 (Case (2)) we must have W3 = which leaves L3 ® / as the only 
constant of motion. However, Case (1) puts no constraint on the z-component of 6«, and 
if the coupling vector is perpendicular to the motion of the plane (b z 7^ 0), both L3 ® I 
and L3 <g) as are constants of motion but of course they are not independent. Therefore, 
including the Hamiltonian there are two constants of motion, independent of which of the 
above cases are satisfied. The first-order integrability is preserved in the presence of a LV 
background b with even arbitrary nonzero components in each direction. 

3.3 Lorentz violation with pure timelike axial- vector coupling (60 7^ 0): 

When a = 0, b = but bo ^ the Hamiltonian Eq. (2.5) and the integral of motion Eq. 
(2.6) reduce to the following in two-dimensions 

\(fi +Pt) + V (p) + Vy{p) (T3L3 - b {<j lPl + o 2 p 2 ) , (3.13) 
(A Q + A ■ a)pi + (B Q + B- a)p 2 + $0 + * • <x 
- l -{(A + A-a) x + (B + B-a) y } , (3.14) 

where Vq = T\/p and V\ = I^/p 3 are assumed. Again, because of the fact that the 
Hamiltonian Eq. (3.13) is not diagonal, we can no longer choose the integral of motion 
diagonal as in the previous case. 

In this case we have 12 equations coming from the coefficients of the second-order 
terms. Three of them are exactly the same with Eq. (3.3) for u = and hence we have the 
same Aq and Bq as in Eq. (3.4). The rest of the determining equations from the coefficients 
of the second-order terms can be written as 

2b A 2 + A 3x = 0, 2b B 1 -B 3y = 0, 2yA 2 V 1 - A lx = , 2x3^ - B 2y = , 
2b A 3 - 2yA 1 V 1 - A 2x = , 26o5 3 + 2x£ 2 Vi + Bii/ = , 2b A 1 - 2b B 2 - A 3y - B 3x = , 
2b A 3 + 2V 1 {xA 2 - yB 2 ) + A ly + B lx = , 2b B 3 + 2V 1 (xA 1 - yB{) - A 2y - B 2x = . 

(3.15) 

There is only a so called trivial solution for the system of equations given in Eq. (3.15), 
which is A = and B = 0. This is shown again by using the Maple software [17]. 

After introducing the already found functions A$, Bq, A, and B into the determining 
equations coming from the coefficients of the first-order terms we have the following 8 
partial differential equations for $0 and * 

$0x = , <S> y = , 
Sis = 2yVi$ 2 , ®iy = -2xVi<$> 2 - b (w + 2$ 3 ) , 
$2x = -2j/Vi$i + 60(^0 + 2$ 3 ) , $2y = 2xVi$i , 
$3x = -y{w y + €o)Vi x + y{w x - rj )V ly - r] Vi - 26 $ 2 , 
<S>3y = x(w y + £o)Vi x - x(w x - r) )V ly + £oV"i + 26 $i . (3.16) 

It is immediately seen that <3?o = const, which can be taken as zero without loss of generality. 
Then we are left with 6 first-order partial differential for 3>. The requirement of the equality 



u-(^O) 
n 2D 

A 2D 
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of the mixed partial derivatives give 3 more equations for $ and their first-order derivatives. 
Now, if we eliminate the first-order derivatives of <& in this system by using the equations 
given in Eq. (3.16) we get a system of algebraic equations for In order to solve this 
system the rank of the coefficient matrix has to be equal to the rank of the extended matrix, 
determined by adding the non homogeneous vector as a column. A simple analysis shows 
that this condition can only be possible if £o = and 770 = 0. In this case we have the 
following solution 

Finally, we get the following 4 partial differential equations from the coefficients of the 
zeroth-order terms 

2b $3x ~ 2Vi(!/$i a . - X$i y ) - {§2xx + <&2yy) = , 26 ($1< / - $2:r) - ($3xx + $3yy) = , 
2& $3j/ ~ 2Vi(y$2x - X§ 2 y) + ($lxx + = , ?/$2 + x(x<S> 2 y ~ V^2x) = , (3.18) 

where A = 0, B = 0, Aq = wo y and Bq = —wo x are used. Together with Eq. (3.17) it is 
possible to show that only $2 = satisfies simultaneously all the 4 differential equations 
in Eq. (3.18). Thus, to sum up we find the following set of solutions in the case if &o / 0: 

A = w y, B = -w x, A = , B = , 
$o = 0, $i = 0, $2 = 0, $ 3 = -w /2. (3.19) 

It is seen from the constant of motion operator in Eq. (3.14) that we have 

X 2D^° ] = w o{ypi - xp 2 ) - 0-3^0/2 

= ~wo{L% + (J3/2) (3.20) 

which is nothing but the total angular momentum J3. It is important to note here that, 
as opposed to the case in 2D where there are no LV terms and hence both L3 and J3 are 
constants of motion, only J3 is a constant of motion in this case ( 60 7^ 0). The fact that 
L3 does not commute with the Hamiltonian is something expected from the form of the 
LV term in Eq. (3.13). However, first-order integrability is still restored. 



4. The analysis of 3D Pauli system in LV background 

In general in the presence of LV terms we have the Hamiltonian and integral of motion given 
in equations Eq. (2.5) and Eq. (2.6) respectively in 3-dimensional Euclidean spaces. In a 
similar fashion as in the two-dimensional case it is also convenient to analyze the problem 
in three separate cases. However, it is advantageous to give the determining equations 
coming from the coefficients of the second- and first-order terms in full generality (that is, 
keeping a, b, and 6q n on zero) before analyzing them in case by case. 
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i) Determining equations coming from the second-order terms: From the diag- 
onal elements it is immediately found that we have the following 6 determining equations 

Aox = , B 0y = , Cq z = , 
Ao y + Bq x = , Aq z + Co x = , Bq z + Co y = , (4-1) 

from which the following solutions are obvious 

A = Pi- a 3 y + a 2 z , 

Bo = 02 + Ot 3 X - OL\Z , 

Co = 03 - a 2 x + aiy . (4.2) 

where ctj and (i = 1, 2, 3) are integration constants. Note that these solutions are inde- 
pendent of the LV background couplings so that they hold for each of the cases discussed 
below. If we compare Eq. (4.2) with Eq. (3.4) we see that we have the following correspon- 
dences: 0i = £o 5 02 = Vo an d «3 = —wo- After introducing the equation (4.2) into the rest 
of the coefficients of the second-order terms and separating the imaginary and real parts of 
the coefficients coming from the off-diagonal elements we are left with an over determined 
system of 18 partial differential equations for Ai, B{, Cj (i = 1, 2, 3). These are, 



2zA 1 V 1 + A 3x + 2b A 2 = , (4.3) 

2yA 1 V 1 + A 2x -2b A 3 = 0, (4.4) 

2xB 2 V 1 + B ly + 2b B 3 = , (4.5) 

2zB 2 V 1 + B a „-2boB 1 = 0, (4.6) 

2xC 3 V 1 + C lg - 2b C 2 = , (4.7) 

2yC 3 V 1 + C 2z + 2b Q C 1 = 0, (4.8) 

2V!(yA 2 + zA 3 ) - A lx = , (4.9) 

2V 1 (xB 1 + zB 3 ) -B 2y = 0, (4.10) 

2V 1 (xC 1 + yC 2 )-C 3z = 0, (4.11) 

2zV 1 {A 2 + B 1 ) + A 3y + B 3x + 2b (B 2 -A 1 ) = 0, (4.12) 

2yV 1 (A 3 + d)+A 2z + C 2x + 26 (^i - C 3 ) = , (4.13) 

2xV 1 {B 3 + C 2 ) + B lz + C ly + 26 (C 3 - B 2 ) = , (4.14) 

2V 1 (xA 1 + yA 2 - zd) - A 3z - C 3x - 2b C 2 = , (4.15) 

2V 1 (xB 1 + yB 2 - zC 2 ) - B 3z - C 3y + 2b d = , (4.16) 

2V 1 (xA 2 - yB 2 - zB 3 ) + A ly + B lx + 2b A 3 = , (4.17) 

2V 1 (xA 1 + zA 3 -yB x )- A 2y - B 2x + 2b B 3 = , (4.18) 

2V!(xA 3 - yC 2 - zC 3 ) + A lz + C lx - 2b A 2 = , (4.19) 

2Vi {yB 3 - xd -zC 3 ) + B 2z + C 2y + 2b B 1 = . (4.20) 
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ii) Determining equations coming from the first-order terms: We set the coeffi- 
cients of the pi (i = 1,2,3) zero at each entry of the commutation relation. After intro- 
ducing equations (4.2) and separating the real and imaginary parts, we have the following 
21 partial differential equations, 12 of which are 



X f 


■0Vi + (Pi- ya 3 + 2y<S> 3 )Vi + $ 2z - a • VC 2 


- b (ai - 


2$i) - 2(6 x C) 2 


= 0, 


(4.21) 


X f 


■0Vi + (Pi + za 2 - 2z$> 2 )V 1 - $ % + a ■ VB 3 


- b (ai - 


2$i) + 2(b x B) 3 


= 0, (4.22) 


yr 


■PVi + ifa- zai + 2z$ 1 )V 1 + $ 3x -a- VA 3 


- b (a 2 - 


2$ 2 ) - 2(b x A) 3 


= 0, 


(4.23) 


yr 


■PVi + (fo + xa 3 - 2x<S> 3 )V 1 -&iz + a- VCi 


- b (a 2 - 


2$ 2 ) + 2(bxC)i 


= o, 


(4.24) 


z f ■ 


PVl + (/% - xa 2 + 2x$ 2 )Vi + *i w - a ■ VB 1 


- b (a 3 - 


2$ 3 ) -2(6x B)i 


= 0, 


(4.25) 


zr ■ 


■0Vi + (03 + yai - 2y$ 1 )V 1 -$ 2x + a- VA 2 


- b (a 3 - 


2$ 3 )+2(fexA) 2 


= o, 


(4.26) 




(ya 2 + za 3 - 2y$> 2 - 2z$ 3 )Vi + $ lx - 


- a ■ VAi 


- 2(6 x A)i = 0, 




(4.27) 




(xai + za 3 - 2x<Z> 1 - 2z$ 3 )V 1 + <& 2y - 


- a ■ VB 2 


-2(bxB) 2 = 0, 




(4.28) 




(xai + ya 2 - 2x$> 1 - 2y<^ 2 )V 1 + $> 3z - 


- a • VC 3 


-2(6xC) 3 = 0, 




(4.29) 



({yA 3x - xA 3y ) + {xA 2z - zA 2x ) + {zA ly - yA lz ) + (C 2 - B 3 ))ti 

- (x(f x A)i + y{f x B)i + z{f x C)i) VI - (a x a)i + 6 V • A , (4.30) 
[{yB 3x - xB 3y ) + (xB 2z - zB 2x ) + (zB ly - yB lz ) + (A 3 - Ci)) Vi 
- (x(f x A) 2 + y(r x B) 2 + z(r x C) 2 ) Vi - (a x a) 2 + b V • B , (4.31) 
{(yC 3x - xC 3y ) + (xC 2z - zC 2x ) + (zC ly - yC lz ) + {B x - A 2 f)V 1 

- (x(f x A) 3 + y(r x B) 3 + z(r x C) 3 ) VI - (a x a) 3 + 6 V • C , (4.32) 

where Vf = dV\/dr, f is the unit displacement vector, a = (a\,a 2 ,a 3 ) and b = (b\,b 2 ,b 3 ) 
are the space parts of the Lorentz violating parameters and b^ respectively, and A = 
(A\, A 2 , A 3 ) (B and C are defined similarly) is introduced in the integral of motion X. 
The rest of the determining equations from the coefficients of the first-order terms can be 
expressed in terms of certain derivative combinations of the determining equations from 
the coefficients of the second-order terms given in Eqs. (4.3) - (4.20). For completeness 
they are symbolically listed below: 



2(Eq(4.9)) x - (Eq(4.17))„ - (Eq(4.19)) 2 = , (4.33) 

2(Eq(4.4)) x - (Eq(4.18)) w + (Eq(4.13)) 2 = , (4.34) 

2(Eq(4.3)) x + (Eq(4.12)) w - (Eq(4.15)) 2 = , (4.35) 

2(Eq(4.5))„ + (Eq(4.17)) x + (Eq(4.14)) 2 = , (4.36) 

2(Eq(4.10)) y + (Eq(4.18)) a - (Eq(4.20)) 2 = , (4.37) 

2(Eq(4.6)) w + (Eq(4.12)) x - (Eq(4.16)) 2 = , (4.38) 

2(Eq(4.7)) 2 + (Eq(4.19)) x + (Eq(4.14))„ = , (4.39) 



$o x = 
<5>0z = 
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2(Eq(4.8)) 2 + (Eq(4.13)) x + (Eq(4.20))„ = , (4.40) 
2(Eq(4.11)) 2 + (Eq(4.15)) a + (Eq(4.16))„ = , (4.41) 

where (...) x for example represents the derivative of the entire equation with respect to x. 

We do not include the expressions for the 8 determining differential equations coming 
from the zeroth-order terms. They are very lengthy and their exact form is not partic- 
ularly illuminating. We simplify the equations by imposing the solutions obtained from 
the second- and first-order terms, and then they are either trivially satisfied or reduce to 
simple constraint equations (like Eq. (3.7) in the 2D case), presented for each case in the 
following discussion. 

In all three cases to be discussed below we have the same integral of motion Eq. (2.6) 
and the functions Aq, Bq and Co are obtained the same as in Eq. (4.2). Therefore, we will 
briefly discuss the solutions for the rest of the coefficients by applying the technique used 
so far. 

4.1 Lorentz violation with vector coupling (eij ^ 0) : 

In this case we have the Hamiltonian Eq. (2.5) with 6 = and bo = 0. It is found by using 
Maple that when bo = we have the following solution for the 18 determining equations 
given in Eq. (4.3) - Eq. (4.20) 

Ai = 0, A 2 = zw, A 3 = -yw, 
B\ = —zw, .E>2 = 0, B3 = xw, 

Ci = yw, C 2 = -xw, C 3 = 0, (4.42) 

where w is an integration constant. 

After introducing Eq. (4.42) into the rest of the determining equations it is immediately 
found from Eq. (4.30) - Eq. (4.32) that 

$0 = (0203 - 0302)2 + ("301 - "103)2/ + (aia 2 - a 2 ai)z , (4.43) 

and then we are left with 9 first-order partial differential equations for 3> = ($1, $2, $3) ■ 
These are the equations Eq. (4.21) - Eq. (4.29) with 6 = and 60 = 0. In order to 
solve this system we express the first-derivatives of and look for the compatibility of the 
mixed partial derivatives. The requirement of the equality of the mixed partial derivatives 
gives another 9 equations for 3> and its first-order derivatives. Now, introducing the first- 
derivatives of which are found from Eq. (4.21) - Eq. (4.29) into this system, we get a 
system of algebraic equations for <I>. In order to have a solution of this algebraic system we 
must have (3 = and w = 0. It is not suprising to require w = since it is the coefficient 
of the a ■ L term in Eq. (2.6) and clearly it cannot commute with the Hamiltonian Eq. 
(2.5) in the presence of a terms. Setting (3 = and w = 0we find 

* = f ( 4 - 44 ) 

Finally, if we introduce all the information gathered from the coefficients of the higher- 
order terms into the determining equations coming from the coefficients of the zeroth-order 
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terms, we find that the following 3 conditions must also be satisfied 

ai(ya 2 + za 3 ) - (a 2 y + a 3 z)a\ = , 
a 2 (xa\ + za 3 ) - {a\x + a 3 z)a 2 = , 

a 3 (xa\ + ya 2 ) - (a\x + a 2 y)a 3 = . (4.45) 

Clearly it is seen that in order to satisfy these 3 conditions we must either have a = or 
a. = 0. Since the former case implies that there are no LV terms and the latter one indicates 
we have no nontrivial first-order integral of motion, we conclude that in the presence of a 
terms there is no first-order integrability unlike the 2D case, if a is perpendicular to the 
plane of the motion (L3 as discussed is constant of motion). It should be noted that 
term can be eliminated with a transformation and in general is not observable physically. 
This is only true for simple one fermion models or non-interacting multi-fermion cases but 
does not hold for a generic multi-fermion theory. 

4.2 Lorentz violation with pure spacelike axial- vector coupling (6j 7^ 0): 

After setting a = and 60 = in Eq. (2.5), we again find the same solution as in Eq. (4.42) 
for A,B, and C since the 18 equations are identical in the absence of bo- However, from 
Eq. (4.30) - Eq. (4.32), in this case we obtain <£ = const, which can be taken as zero 
without loss of generality. Then, we proceed in a similar fashion as in the previous case 
and reach the same conclusion for the system of 9 differential equations given in Eq. (4.21) 
- Eq. (4.29), with a = and bo = 0. That is, $ = a/2 is still the solution together with the 
requirement that (3 = and w = 0. After introducing all the information gathered from 
the coefficients of the higher-order terms into the determining equations coming from the 
coefficients of the zeroth-order terms, we find, as in the previous case, that the following 
condition must also be satisfied 

b ■ a = . (4.46) 

Again, we can conclude that in general in the presence of b terms there is no first- 
order integrability. However, we may find some integral of motions in some special cases 
(e.g. choosing b± = 0, b 2 = and a± = 0, a 2 = 0). In general there is a constant of 
motion = a ■ J if b is a vector perpendicular to a ( here J represents the total 

angular momentum operator). Of course, not all components of en are independent due to 
Eq. (4.46). For example, if we assume the i th component of b to be nonzero, then one can 
alternatively write the constant of motion as 

4T ] = r[(«x(*x J)\ ■ (4.47) 

where a is otherwise arbitrary. From here one could think that there are two constants 
of motion after eliminating one of the components of a but this is not true. Once the 
orientation of b is fixed, there is a unique a so that we have one additional constant of 
motion, which could be chosen as J3. One noticeable difference from the 2D case is that the 
total angular momentum J3 is now conserved, but this is not enough to make the system 
first-order integrable. 
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4.3 Lorentz violation with pure timelike axial- vector coupling (bo ^ 0): 

In this case, the set of 18 equations in Eqs. (4.3) - (4.20) differs from the previous 2 cases 
since &o ^ 0. The only solution of this system is the so called trivial solution which is 
A = 0, B = and C = 0. We find also that $o = const, which can be taken as zero 
without loss of generality, and proceed in the same way as in the previous cases. We again 
find the same <& as in Eq. (4.44) together with the only requirement (3 = 0. However, this 
time the determining equations coming from the coefficients of the zeroth-order terms are 
identically satisfied if we introduce all the information gathered from the coefficients of the 
higher-order terms. The form of X^^ ^ is the same as in the previous case but without 
the constraint Eq. (4.46). There exist however two additional constants of motion since 
the third one is related to the first two through an SU(2) algebra. 

Thus, we conclude that the Hamiltonian Eq. (2.5) with a = and b = is first-order 
integrable even though it has LV terms (6o ^ 0). It is easily seen that the components 
of J commute with H, hence the system is integrable. We may choose {Jj, J 2 , H} as a 
commuting set of operators. It is important to note that in the commuting set we choose 
J 2 but there is no way to get it in since J 2 is a second-order operator. However, 

first-order integrals of motion Jj generate J 2 . 

5. Perturbative Solutions 

The whole idea in this study is to seek, with a systematic method, the first-order inte- 
grability/superintegrability of the Pauli system in (axial-)vectorial LV background. The 
ultimate goal is to achieve at least a quasi-exact solution. This could be important for the 
6o 7^ case since the bound on bo for electron is much weaker than both a and b, and if 
the effective size of the system comparable with 6q 1 , the perturbative approach fails. So, it 
would be enough to mainly concentrate on the bo ^ case and even perform a second-order 
analysis. However, our experience is that it would be much harder to solve the equations 
of overdetermined systems analytically, especially if one keeps the spin-orbit term. 

Based on our analysis we don't have first-order superintegrability in either 2D or 3D. 
There are first-order integrable cases in 2D under special arrangements but this is only true 
in 3D for the bo ^ case, where exact solutions could be of interest. Of course, neither 
integrability nor superintegrability does guarantee the separability and existence of exact 
solution but either is an important step in finding one. Recently, a perturbative approach 
has been carried out for a similar Hamiltonian system in the same LV background within 
the QED version of Standard Model Extension [18, 19]. For completeness, we like to briefly 
summarize their results. Note that the analysis in Ref. [18] includes an external magnetic 
field but not the spin-orbit term. A more extensive perturbative analysis for all types of 
Lorentz violating couplings is found in [20]. 

For the case with a^, as discussed in Section 4.1, using a field redefinion by a factor of 
exp(i a^x 11 ), the energy is shifted by a constant ao and wavefunction has an overall phase 
factor exp(ia ■ r), compared to a hydrogen atom without LV background. 

For the case 6^0, the perturbative approach is applied [18] and energy shifts are 
obtained, by an amount proportional to b^m s (to first order). Here m s is the magnetic spin 
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quantum number. The wave function is unchanged. The energy correction is proportional 
to b-$rrij/(l + 1) if the spin-orbit correction is included. Here rrij and I are the total 
magnetic angular momentum and the orbital angular momentum numbers, respectively. 
These corrections are quite small due to experimental bounds on |6|. 

The case with bo / is indeed interesting. It is shown that using a unitarity transfor- 
mation one could relate a real hydrogen with all corrections (spin-orbit, relativistic, and 
the so-called Darwin term) to the one with additional LV background of the form b^a • p. 
The unitarity transformation is U = 1 — i bo(l + S/r)<r • r, up to first order in bo, so that 
^rcamy = ^ H TealHy U and ^( b o^°) = . Here S is a constant appropriately chosen 
to obtain the correction terms. It is then straightforward to solve the system i? rea iHy by 
applying perturbation theory and to transform everything back to the original system. As 
usual there is no contribution to the energy spectrum of the hydrogen atom from the bo 
term, up to first order in bo- However, the wave function develops a bo dependent part. 
See the Ref. [18] for the details. 

6. Conclusion 

In this study we have analyzed in a systematic way the integrability of a Pauli system in 
a LV background both in two- and three-dimensions. We consider the Hamiltonian which 
originates from the non-relativistic limit of the QED part of the so-called Standard Model 
Extension. We kept only two types of LV couplings to the electron, the vector type of 
coupling represented by a M and the axial- vector type of coupling by b^. Since ao enters as 
a constant term in the Hamiltonian, we consider in both 2D and 3D the following three 
separate cases; LV from a / 0, LV from 6^0, and LV from &o 0- The last case is 
especially important since perturbative approaches could fail under certain circumstances, 
due to weak experimental bounds on bo, and exact solutions might be required. Note that 
in our study we kept the spin-orbit correction to the hydrogen atom, partially because the 
correction from the LV term with bo / could be comparable with the spin-orbit term, as 
shown in Section 2. 

After writing the most general constant of motion operator X in first-order deriva- 
tives, we obtained sets of overdetermined differential equations from the vanishing of the 
commutator of X with the Hamiltonian. In 2D, first-order integrability is possible if a is 
perpendicular to the plane of the motion, when the LV is due too^O only. The constant 
of motion is one of the components of the angular momentum operator L (£3). Our con- 
clusion is the same for b 7^ case, with the difference that the vector b can have arbitrary 
components in each direction. The last case in 2D is 60 7^ and first-order integrability is 
restored like in the other cases with one difference: the z-component of the total angular 
momentum operator J3 is conserved if the motion is in xy plane. 

Based on our 2D results, we extended our discussion into a more realistic 3D picture 
by following the same procedure. The first-order integrability disappears for the o / 
case in 3D. For 6^0 case, even though J3 becomes a constant of motion, it is not enough 
to make the system integrable (since we need another one, unlike 2D case). The first-order 
integrability in 3D is retained only for the bo 7^ case where we can find two additional 
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constant of motion operators like J\ and J2. We believe that studying the 60 / case 
up to second-order is worthwhile in search of superintegrability. This is in fact the case 
where exact solutions could be physically relevant. One caveat is that, at the second order, 
solving the overdetermined partial differential equations becomes much more challenging. 
For completeness, we also summarized the results of a recent paper [18] where unitary 
transformations and perturbative approaches are employed for similar systems. 
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